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ABSTRACT 

To interpret H 2 quasar absorption line observations in Damped Lya clouds (DLAs), 
we model the H2 spatial distribution within a DLA. Based on numerical simulations 
of disk structures with parameters similar to those derived for such absorbers, we 
calculate the H 2 distribution as a function of ultraviolet background (UVB) intensity 
and dust-to-gas ratio. For typical values of these two quantities we find that the area in 
which the H 2 fraction exceeds 10~ 6 (typical observational detection limit) only covers 
< 10% of the disk surface, i.e. H 2 has a very inhomogeneous, clumpy distribution even 
at these low abundance levels. This explains the relative paucity of H 2 detections in 
DLAs. We also show the dependence of the covering fraction of H 2 on dust-to-gas 
ratio and UVB intensity and we comment on the physics governing the H 2 chemical 
network at high redshift. 

Key words: ISM: molecules — galaxies: evolution — galaxies: high-redshift — cos- 
mology galaxies — quasar absorption lines 



1 INTRODUCTION 

In the recent years, evidence has been found for the exis- 
tence of heavy elements in damped Lya clouds (DLAs), i.e. 
quasar absorption line systems whose neutral hydrogen col- 
umn density is larger than ~ 1 — 2 x 10 20 cm -2 (e.g. Pettini 
et al. 2001; Prochaska & Wolfe 2002). The evolution of metal 
abundance in DLAs can trace the chemical enrichment his- 
tory of present galaxies. Based on this, and on other clues, 
DLAs have been suggested to be the progenitors of nearby 
galaxies; the similar values of the baryonic mass density in 
DLAs around redshift z ~ 2 and the stellar mass density at 
z ~ has further supported this idea (Lanzetta, Wolfe, & 
Turnshek 1995). 

It is natural to consider that the DLAs contain a certain 
amount of dust, because dust-to-gas ratio and metallicity are 
correlated even for metal-poor galaxies (Schmidt & Boiler 
1993; Lisenfeld & Ferrara 1998). Indeed, Fall, Pei, & McMa- 
hon (1989) have suggested that the reddening of background 
quasars indicates typical dust-to-gas ratios of ~ 1/20-1/4 of 
the Milky Way (see also Zuo et al. 1997). The depletion of 
heavy elements also supports the dust content in DLAs (e.g. 
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Vladilo 2002). The existence of dust implies the possibility 
that the formation of hydrogen molecules (H 2 ) is enhanced 
because of the H 2 grain surface reaction. Hirashita & Fer- 
rara (2002) have recently shown that even in metal poor 
(~ 0.01 Zq) galaxies, dust grains can drastically accelerate 
the formation rate of H2. They also argued that the enhance- 
ment of molecular abundance results in an enhancement of 
the star formation activity in the early evolutionary stages 
of galaxy evolution, because stars only form in molecular 
clouds. The important role of dust on the enhancement of 
the H 2 abundance is also suggested by observations of DLAs 
(Ge, Bechtold, & Kulkarni 2001) and in the local Universe, 
e.g., in Galactic halo clouds (e.g. Richter et al. 2001) and 
in the Magellanic Clouds (Richter 2000; Tumlinson et al. 
2002). 

Although the H 2 fraction, /h 2 = 2A(H 2 )/[2A(H 2 ) + 
AT (Hi)], where N(X) indicates the column density of the 
species X, is largely enhanced for some DLAs, stringent up- 
per limits are laid on a significant fraction of DLAs in the 
range ~ 10" 7 -10" 5 (Petitjean, Srianand, & Ledoux 2000). 
This can be interpreted as due to a low formation rate of H 2 
in dust-poor environments relative to the Milky Way (Lev- 
shakov et al. 2002; Liszt 2002) and high H 2 dissociation rate 
by strong ultraviolet background (UVB) radiation (e.g., Pe- 
titjean et al. 2000). However, we should keep in mind that 
such upper limits do not exclude the existence of molecular- 
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rich clouds in these systems, because molecular clouds may 
have a very low volume filling factor. If the covering fraction 
of molecular-rich regions on a galactic surface is extremely 
small, it is natural that H2 is hardly detected in DLAs. Thus, 
the argument on the H2 abundance in DLAs is strongly de- 
pendent on the geometry of H2 distribution within those 
systems. 

In order to get a better understanding of the spatial 
distribution of H2, which can then be used to interpret ob- 
servations, we present here a study of such problem based 
on high-resolution numerical simulations. This allows us to 
tackle the problem of the H2 formation/destruction and dis- 
tribution in DLAs in a realistic way. We calculate the spatial 
structure of H2 distribution in a galactic disk under various 
conditions by varying the UVB intensity and dust-to-gas 
ratios. As the underlying gas density and temperature dis- 
tribution is, to a first approximation, independent of the H2 
properties (for example, it does not contribute to cooling at 
the metallicity level typical for DLAs), this approach offers 
the opportunity to explore the parameter space with the 
required accuracy. 

Throughout this Letter, we assume a flat ACDM cos- 
mology 1 (Mo & White 2002). We first describe the simula- 
tion that we used to derive the density and temperature 
maps (Section The maps are presented in Section [3] 
where the H2 distribution maps are also shown. Based on 
these results, we discuss our interpretation of current obser- 
vations of H2 in DLAs (Section 2J. 



2 NUMERICAL SIMULATION OF DLA DISKS 

It is still unclear whether DLAs are large protogalactic disks 
(Prochaska & Wolfe 1998), protogalactic clumps (Haehnelt, 
Steinmetz, & Rauch 1998; Ledoux et al. 1998), or a mix- 
ture of various populations (Cen et al. 2002). Here, we as- 
sume that DLAs are large protogalactic disks. We use a two- 
dimensional hydrodynamical simulation assuming a disk-like 
geometry (for the method, see Wada & Norman 2001). How- 
ever, the statistical properties of density and temperature 
are determined by non-linear hydrodynamical effects (Wada 
& Norman 2001), and we expect that the results of this Let- 
ter remain valid if other scenarios for DLA formation are 
adopted. 

We have run the hydrodynamical calculation code de- 
scribed by Wada & Norman (2001) to obtain density and 
temperature spatial distributions. The parameters for the 
simulation are set as follows. The velocity dispersion derived 
from the line width is roughly 100 km s" 1 (ranging from 40 
to 300 km s _1 ; e.g., Prochaska & Wolfe 1998). Therefore, we 
fix the circular velocity as v c = 100 km s _1 . We assume a 
formation redshift z v j r = 3. By using the spherical collapse 
model for galaxy formation (see equations 2 and 7 of Hi- 
rashita & Ferrara 2002), those values yield a virial mass of 
M v i r = 8.0 x 10 10 Mq and a radius of the dark matter halo of 
r v i r = 34 kpc. If the baryon fraction is assumed to be equal 
to Qb /Om , the gas mass contained in the galactic disk within 
the halo is estimated to be Mdisk = 1.1 X 10 10 /di s k Mq, where 

1 We use the following values of the cosmological parameters: 
tt M = 0.3, n A = 0.7, H = lOO/i km s" 1 Mpc" 1 = 70 km s" 1 
Mpc" 1 , and C b = 0.02h" 2 . 



/disk is the gas fraction contained in the disk. According to 
Navarro & Steinmetz (2000), we estimate /disk — 0.2 for our 
cosmological parameters. We do not include star formation 
and stellar feedback. The UV radiation for the photoelec- 
tric heating is assumed to be 1/100 of the local Galactic 
value. The following conclusions, however, are not affected 
by the details of the local radiation field, because the den- 
sity and temperature distributions of the dense gas where 
H2 molecules are mainly formed are not very sensitive to 
the UV intensity. We have confirmed that the dense regions 
survive even in the strongest UV intensity treated in this 
Letter. 

To solve our problem it is necessary to follow the evo- 
lution of the system from galactic scales down to the small 
regions where H2 forms. We use a resolution of 2048 x 2048 
cells with a fixed cell size of 0.49 pc, so as to reach a global 
1 kpc scale simulation. The typical radius of the gaseous 
disk can be estimated to be 0.18r v i r = 6.1 kpc (e.g., Ferrara, 
Pettini & Shchekinov 2000). With this size, the column den- 
sity of hydrogen nuclei in the vertical direction of the disk is 
~ 1 x 10 22 /disk cm -2 . We simulated the central 450 pc radius 
of the disk, in which we assume that a fraction (0.45/6. 1) 2 
of A/disk is contained. Then, we simulate an exponential disk 
with a scale length of 100 pc, to check the effects of radial 
density profiles. From a series of tests, we have been able to 
assess that the following results are not affected significantly 
by the specific profile assumption. 

The whole timescale of the simulation is 60 Myr, when 
a quasi-stationary density distribution function is achieved. 
This timescale is larger than the rotation time in the sim- 
ulated region: 2tt(450 pc)/100 km s" 1 = 30 Myr. For the 
cooling function, we assume a metallicity of 0.1 Zq, appro- 
priate for DLAs. 

In Fig. we show the density and temperature distri- 
butions calculated by the simulation. A major part (86%) 
of the disk is covered with regions with nu < 10 cm -3 (n,H 
is the number density of hydrogen nuclei). This density is 
roughly consistent with that derived observationally by Silva 
& Viegas (2002) . The gas temperature T is higher than 500 
K in a large part of the disk (62%), and indeed Chengalur 
& Kanekar (2000) observationally derived T > 500 K for a 
large part of their sample. There are also a few DLAs de- 
tected with T ~ 100 K (e.g. Chengalur & Kanekar 2000). 
For more detailed statistical study, a complete set of simu- 
lations that cover a range of mass, rotational velocity, etc. 
is necessary. The above qualitative agreement, however, as- 
sures a good basis on which we discuss the properties of 
DLAs. Those density and temperature maps are used to es- 
timate the spatial distribution of molecular fraction in the 
following. 



3 H 2 SPATIAL DISTRIBUTION 

3.1 Estimate of molecular fraction 

We estimate the molecular fraction (/h 2 ) at each grid point. 
Since we are interested in the typical metallicity of DLAs 
(~ 1/10 Zq; Vladilo 2002), the gas phase formation of H2 
through the formation of H~ is negligible compared with 
the formation on dust (Liszt 2002). On the other hand, the 
dominant mechanism for the destruction of H2 is the pho- 
todissociation under a UVB intensity. 
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Figure 1. Spatial distribution of density (left) and temperature (right) obtained from the simulation. The color bars show the levels of 
log riH and logT, respectively. 



For simplicity, we assume equilibrium between H2 for- 
mation and destruction. In order to examine if this assump- 
tion is justified, we first estimate the formation timescale 
of H2 as fH 2 n a/Ri ~ 10 4 yr, where we adopt the reaction 
rate Ri ~ 4.6 x 10~ 16 cm -3 s _1 (defined in equation 0, 
/h 2 ~ 10 -6 (about the experimental detection limit), and 
riH ~ 100 cm -3 (typical density for molecule-forming re- 
gions) . In this Letter, we are interested in structures down to 
0.5 pc (grid size), and the hydrodynamical timescale in such 
structures can be estimated as 0.5 pc/10 km s" 1 ~ 5 x 10 4 
yr. We are particularly interested in the cold clouds where 
H2 can form, and in such a region, the typical hydrodynam- 
ical timescale is > 10 yr because of their low sound speed. 
Thus, the chemical equilibrium can be reasonably assumed 
for /h 2 ;$ 10 -6 and nu > 100 cm -3 . Since such an equi- 
librium becomes a bad approximation for high fu 2 , some of 
the molecular rich clouds with fu 2 > 10~ 4 will disappear in 
the course of hydrodynamical evolution. In a diffuse gas with 
nu < 10 cm -3 , the chemical equilibrium may not be realised 
because of a slow reaction rate. This suggests that H2 does 
not form so efficiently as estimated from the equilibrium 
assumption in a diffuse medium. Moreover, such a diffuse 
medium has a temperature > 10 3 K. With such a high tem- 
perature, hydrogen atoms may not stick to dust efficiently. 
Therefore, /h 2 in diffuse region may be overestimated in 
this Letter. Fortunately this strengthen our conclusion on 
the lack of H2 in a diffuse medium. A more consistent treat- 
ment of H2 reaction and hydrodynamics will be tackled in 
the future. 

We adopt the formation rate of H2 per unit volume and 
time, Ri, by Hollenbach & McKee (1979) (see also Hirashita, 
Hunt, & Ferrara 2002): 



Ri — 0.5ti,h(1 — fu 2 )nd^a vSd(T) 
~ 4.6 x 10 _16 5 d (T) 
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a grain (assumed to be spherical with a radius of 0.1 nm), v 
is the mean thermal speed of hydrogen, T> is the dust-to-gas 
mass ratio, 5 is the grain material density (assumed to be 
3 g cm -3 ), and Sd(T) is the sticking coefficient of hydrogen 
atoms onto dust. In equation we have substituted the 
dust number density rid with the dust-to-gas ratio T> by 
using 

rid— a 3 5 = numuV , 

where mu is the mass of a hydrogen atom. The sticking 
coefficient is given by (Omukai 2000) 

S d (T) = [l + 0.4(T + T d ) ' 5 +2 x 10" 3 T + 8 x lO^T 2 ]" 1 



(2) 



x[l + exp(7.5 x 10 2 (l/75- 1/T d ))]" 



(3) 



where Td is the dust temperature. In this Letter, T d is as- 
sumed to be 20 K (a typical temperature under the local 
interstellar radiation field), but the following result is insen- 
sitive to this value as long as Td < 70 K. 

The photodissociation rate in units of cm -3 s -1 , 7?2, is 
estimated by (Abel et al. 1997) 

R 2 = (4tt) 1.1 x 10" 13 n H2 J21 Shield [iV(H 2 ), JV(dust)] , (4) 

where iih 2 (cm" 3 ) is the number density of H 2 , J21 (10 -21 
erg s _1 cm -2 Hz" 1 str" 1 ) is the UVB intensity at the Ly- 
man limit wavelength (912 A) averaged for all the solid an- 
gle, S s hield(iV(H2), TV(dust)) is the correction factor of the 
reaction rate for H2 self-shielding and dust extinction. We 
adopt the correction for the H2 self-shielding by Draine & 
Bertoldi (1996) (see also Hirashita & Ferrara 2002). Then, 
we estimate 5*shieid as 



S, 



shield = mm 



f nnfn 2 H V 
V 10 14 cm" 2 J 



(5) 



where rid is the number density of grains, a is the radius of 



where H is the typical thickness of the disk which is assumed 
to be 100 pc in this Letter. 

In fact, in some DLAs with H2 detection, the H 2 excita- 
tion is consistent with the radiation field comparable to the 
local Galactic radiation field (Ge & Bechtold 1997; Ledoux, 
Srianand, & Petitjean 2002; Petitjean, Srianand, & Ledoux 
2002). In this case, the photo-dissociation rate can be much 
higher. However, the H2-detected DLAs could be biased to 
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star- forming molecular regions, and it is still unknown if 
DLAs in general are exposed to such a strong radiation field. 
Thus, we concentrate on the UVB, which is common for all 
the DLAs. The equilibrium condition Ri = R2 with tem- 
perature and density at each grid point gives /h 2 for each 
point. Thus, we obtain the spatial distribution of /h 2 under 
a certain set of (T>, J21). 



3.2 Molecular fraction maps 

The metallicity level of DLAs (~ 0.1 Zq) implies that the 
dust-to-gas ratio of DLAs is typically 10% of the Galactic 
(Milky Way) value. In this Letter, we assume the Galactic 
dust-to-gas ratio to be 0.01. Therefore, we examine the dust- 
to-gas ratio around 0.001. The UVB intensity is considered 
to be J21 ~ 0.3-1 around z ~ 3 (Giallongo et al. 1996; 
Cooke, Espey, & Carswell 1997; Scott et al. 2000; Bianchi, 
Cristiani, & Kim 2001). Since the UVB intensity is likely 
to be lower at z < 1 (e.g. Scott et al. 2002), we examine 
J21 ~ 0.01-1 in this Letter. 

In Fig. |5] we show the molecular fraction (/h 2 ) for (a) 
T> = 0.001 (10% of the Galactic dust-to-gas ratio and J21 = 
0.3, (b) V = 0.001 and J21 = 1, and (c) V = 0.03 and 
J21 = 0.1. The contour levels are for /h 2 = 10" 6 , 10" 5 , 
and 10 -4 . In order to see the fine structure of the molecular 
distribution, we also zoom on the region around (700 pc, 
700 pc) in Fig. |3] We find that molecular-rich regions are 
distributed very inhomogeneously and concentrate in small 
(~ 10 pc) clumps. 

We define the covering fraction C(> /) as the fraction 
of the area where the molecular fraction is larger than / in 
the central 450 pc (~ radius of the simulated area) region. 
The disk is seen face-on and projection effects of disk inclina- 
tion is not considered in this paper. Since the disk thickness 
is much smaller than the size of the galaxy, the effect of 
the inclination on the surface covering fraction is expected 
to be small. Such a projection effect can be quantitatively 
investigated in a future three-dimensional simulation. 

Since the typical detection limit for /h 2 is 10~ 6 , we 
are particularly interested in C(> 10~ 6 ). Observationally, 
C(> 10~ 6 ) can be interpreted as the probability of detecting 
a gas with /h 2 > 10 -6 under a condition that such a disk 
as simulated in this Letter hits the line of sight to a distant 
quasar. If we assume that the simulated disk is typical for 
DLAs, C(> 10~ 6 ) approximates the probability of detecting 
H 2 in a given DLA with a level of /h 2 > 10" 6 . In Tabled we 
show C(> 10 -6 ) as a function of T> and J21. The covering 
fraction is very sensitive to both parameters in the range of 
interest. For another observational thresholds, a power-law 
scaling C(> /) cx f~ a is applicable in 10 -7 < / < 10~ 4 
with a — 0.11, 0.13, and 0.10, for the same parameter sets 
as Figs. [5^,, b, and c, respectively. 



4 SUMMARY AND DISCUSSION 

Our results indicate that the fraction of the area with de- 
tectable molecular fraction (typically /h 2 > 10 -6 ) is very 
small for small dust-to-gas ratio {T> < 0.001, i.e. < 10% of 
the Galactic dust-to-gas ratio) and in a strong UVB radia- 
tion ( J21 > 0.1). Thus, it is rare that a molecular-rich region 
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Figure 2. Spatial distribution of molecular fraction (/h 2 ) for (a) 
T> = 0.001 (10% of the Galactic dust-to-gas ratio) and J21 = 0.1, 
(b) V = 0.001 and J 2 i = 0.3, and (c) V = 0.003 and J 2 i = 0.1. 
The contour levels arc for / Ha = 10~ 6 , 10~ 5 , and 10~ 4 (blue, 
green, and red, respectively). 



hits the line of sight to a distant quasar. This means that 
DLAs with a detectable molecular fraction are rare objects. 

The lack of H2 detection from DLAs does not necessar- 
ily mean the lack of gas in molecular form. Under a strong 
UVB, molecules are predominantly confined in small areas. 
However, in such areas, the molecular fraction is as high as 
> 10 -3 . This molecular fraction is large enough to provide 
the minimum radiative cooling necessary to ignite the star 
formation process. 
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Figure 3. Zooming on the region around (700 pc, 700 pc) in Fig. 
l2*k . The color bars show the levels of log /jj 2 • 



Table 1. Covering fraction of the area with /jj 2 > 10 6 (above 
the typical observational detection limit) in % for various T> and 
J21 



J21 


V 0.0001 


0.0003 


0.001 


0.003 


0.01 


0.01 


3.8 


15 


34 


57 


85 


0.03 


1.7 


8.5 


25 


46 


72 


0.1 


0.75 


3.7 


16 


34 


60 


0.3 


0.43 


1.8 


9.5 


25 


48 


1 


0.26 


0.81 


4.9 


16 


37 



Table Q shows that the covering fraction of molecular- 
rich regions sensitively changes as the dust-to-gas ratio 
and/or the UVB intensity vary. If the dust-to-gas ratio is 
10% of the Galactic value and the UVB is as strong as ex- 
pected at high redshift (J21 > 0.1), C(> 10 -6 ) < 10%; on 
the contrary C(> 10 -6 ) > 20% if the dust-to-gas ratio is 
> 30% of the Galactic value. Such sensitive dependence on 
T> may explain the observational correlation between molec- 
ular fraction and the dust abundance for DLAs (Ge et al. 
2001), similar to the local Universe. This underlines the im- 
portance of dust for the H2 chemical network at high red- 
shift. 

We also found that a detectable amount of H2 is lo- 
calised in regions with riH > 100 cm -3 and T < 100 K. 
Those ranges of density and temperature are consistent with 
those observationally derived for Eb-detected DLAs (e.g., 
Ledoux et al. 2002). Although most of the DLAs may con- 
tain H2, molecular clouds are difficult to detect in absorption 
if they are embedded in a dust-poor and/or UV- intense en- 
vironment, because a large fraction of the area is covered 
by a diffuse, molecule-poor medium. Indeed, most of DLAs 
are observationally suggested to arise selectively in diffuse 
neutral gas not associating with a detectable amount of H2 
(e.g. Petitjean et al. 2000). Therefore, the observed lack of 
molecules from absorption measurements does not neces- 
sarily indicates the lack of molecular clouds, but probably 
reflects the small size of molecular-rich regions. 
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